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ABSTRACT 

The enrichment of Pb in the Galaxy is followed in the framework of a detailed 
model of Galactic chemical evolution that already proved adequate to reproduce 
the chemical enrichment of O and of the elements from Ba to Eu. The stellar 
yields are computed through nucleosynthesis calculations in the Asymptotic Gi- 
ant Branch (AGB) phase of low- and intermediate-mass stars, covering a wide 
range of metallicities. The physical parameters of the stellar structure were de- 
rived from full stellar evolutionary models previously computed. We show that 
low-mass AGB stars are the main producers of Pb in the Galaxy, with a complex 
dependence on metallicity and a maximum efficiency at [Fe/H] ~ —1. Our calcu- 
lations succeed in reproducing the abundances of Pb isotopes in the solar system: 
the role attributed by the classical analysis of the s process to the strong com- 
ponent, in order to explain more than 50% of solar ^°^Pb, is actually played by 
the high production of Pb in low-mass and low-metallicity AGB stars. We then 
follow the Galactic chemical evolution of Pb isotopes and give our expectations 
on the s-process contribution to each of them at the epoch of the solar system 
formation. Finally, we present new spectroscopic estimates of Pb abundance on a 
sample of field stars and compare them, together with a few other determinations 
available, with the predicted trend of [Pb/Fe] in the Galaxy. 
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Subject headings: nucleosynthesis - stars: abundances, AGB and post-AGB - 
Galaxy: evolution, abundances 

1. Introduction 

Neutron captures are the main way for the building up of the nuclei heavier than iron. 
According to the seminal work by Burbidge et al. (1957), in order to reconstruct the Galactic 
evolutionary history of heavy elements one has to consider two major mechanisms of neutron 
addition: the s-process (slow neutron-captures) and the r-proccss (rapid neutron-captures). 
For the s-process, the abundance distribution in the solar system was early recognized to 
derive from a non unique site. Indeed, the phenomenological approach of s-process studies, 
or classical analysis, was shown to require three components to account for the abundances 
of all s-nuclei from Fe to Bi: the weak, the main, and the strong component (Seeger, Fowler, 
& Clayton 1965; Clayton & Rassbach 1967; Kappeler et al. 1982; Kappeler, Beer, & Wis- 
shak 1989). 

The weak component, responsible for the s-process nuchdes up to A ~ 90, was ascribed 
to neutron captures occurring in advanced evolutionary phases of massive stars, through the 
activation of the ^^Ne(Q;,n)^^Mg reaction (Lamb et al. 1977; Prantzos et al. 1990; Raiteri 
et al. 1993 and references therein). The main component, feeding the heavier s-process 
nuclides up to A = 208, is known to originate in AGB stars of Galactic disk metallicity 
during recurrent thermal pulses, hereafter TP. Many theoretical and observational works 
(see Busso, Gallino, & Wasserburg 1999 for a review) converge on the idea that Low-Mass 
Stars (hereafter LMS, M < 4 Mq) play the dominant role in s-processing, and that neutrons 
are released in radiative conditions in the interpulse phase via the ^^C(Q;,n)^^0 reaction 
(Straniero et al. 1995). A tiny ^^C pocket is thought to develop in the top layers of the 
He- and C-rich zone (called He intershell) as a consequence of the penetration of a small 
amount of protons from the envelope, to burn locally and to produce neutrons before the 
development of the next TP. During advanced TPs a second small neutron burst is released 
by the marginal activation of the ^^Ne(Q;,n)^^Mg reaction. Finally, the strong component 
was advanced on theoretical grounds by the classical analysis of the s-proccss (Clayton & 
Rassbach 1967), in order to reproduce more than 50% of ^°^Pb, the most abundant Pb isotope 
in the solar system, although the astrophysical site for this process remained unknown. 

In this work we study the Galactic enrichment of Pb in the light of AGB nucleosynthesis 
results presented by Straniero et al. (1997) and Gallino et al. (1998). One of the main goals 
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of these works was the development of a new set of model calculations for AGB stars with 
different masses and metallicities, to provide stronger evolutionary supports to nucleosynthe- 
sis calculations. The results demonstrated in detail the dependence of the s-process yields on 
stellar metallicity, and introduced the idea that a separate astrophysical site for the strong 
component might be unnecessary (Gallino et al. 1998), because a large production of ^osp^ 
derives from AGB stars of low metallicity. This expectation must however be substantiated 
with a detailed model for the Galactic enrichment, to verify whether the production of Pb 
isotopes is really adequate. 

The paper is organized as follows: in § 2 we discuss the sensitivity of Pb production (in 
particular ^°^Pb) to metallicity and to different assumptions on the ^^C pocket in low-mass 
AGB stars. We also analyze the role of intermediate-mass AGB stars in the Pb production at 
different metallicities. In § 3 we present spectroscopic estimates by synthetic spectra of the Pb 
abundance on a limited sample of field metal-poor stars using the Pb line at 3683.48 A and 
outline the difficulties and the aims of such an analysis. In § 4 we briefly introduce the 
Galactic Chemical Evolution (hereafter GCE) model adopted, present our predictions for 
the s-components of the Pb isotopes in the solar system, and compute the [Pb/Fe][] trend 
vs. metallicity in the three Galactic zones in which the GCE model is organized (halo, 
thick disk and thin disk). Then we compare our model predictions with the above new 
spectroscopic determinations and with a few other published data, outlining the need of a 
dedicate observational effort. Finally, in § 5 we summarize the main conclusions and point 
out a few aspects deserving further analysis. 



2. Pb production by s-process nucleosynthesis 

We compute stellar yields for s-processing of Pb in AGB stars as a result of post- 
process calculations on the basis of stellar evolutionary models obtained with FRANEC 
(Frascati Raphson-Newton Evolutionary Code; see Chieffi & Straniero 1989) and presented 
in Straniero et al. (1997) and Gallino et al. (1998). Additional models have been discussed 
by Vaglio et al. (1999) and Straniero et al. (2000). The detailed evolutionary computations 
were performed at [Fe/H] = 0.0, —0.3 and —1.3 for the AGB stars of mass 1.5 and 3 Mq. 
For the 5 Mq case the metallicities were [Fe/H] = 0.0 and —1.3, while for the 7 Mq star only 
models of solar metallicity were calculated. The s-process results has then been obtained 
with post-process runs span over the metallicity interval from Z = Zq down to Z = Z0/4OOO 
for 1.5, 3, 5 and 7 Mq. The results of the FRANEC stellar structures are extrapolated to 



^We adopt the usual spectroscopic notation [Pb/Fe] — log (Pb/Fe) — log (Pb/Fe)© 
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TABLE 1 
Cumulative mass (in Mq) 

DREDGED-UP FROM THE HE INTERSHELL 
TO THE ENVELOPE OF AGB STARS 



[Fe/H] 1.5 Mq 3 Mq 5 Mq 7 Mq 







1.3 


0.07 


0.16 


0.20 


0.20 


0.5 


0.04 


0.09 


0.20 


0.20 


0.0 


0.02 


0.05 


0.10 


0.10 



the whole metaUicity range adopted for the s-process calculations (more details are available 
in Straniero et al. 1997, 2000, and GaUino et al. 1998). Mass loss was taken into account by 
adopting the parameterization of Reimers (1975). The values adopted for this work are rj — 
0.3 for 1.5 Mq , and r) = 1.5 for 3 Mq (Straniero et al. 1997). Calculations for 5 Mq and 7 
Mq are instead performed with rj = 10 and rj = 3, respectively (Straniero et al. 2000). The 
models self-consistcntly reproduce the so-called third dredge up episode (TDU), defined as 
the penetration of the convcctivc envelope below the H-Hc discontinuity. TDU was found 
to occur after a limited number of TPs and to cease when the envelope mass decreases by 
mass loss below a limiting mass of ~ 0.5 Mq. 

The cumulative mass for the He intershell dredged up into the envelope of AGB models 
of 1.5, 3, 5, and 7 Mq according to FRANEC models and eventually ejected into the inter- 
stellar medium by stellar winds is reported in Table 1. Different treatments for convection 
with respect to the one adopted in FRANEC, where the convective/radiative borders are 
followed according to the Schwarzschild criterion, may provide larger TDU efficiencies, as e.g. 
discussed by Frost & Lattanzio (1996), Herwig et al. 1997; Langer et al. 1999; Herwig 2000). 
Moreover, the criterion adopted for the mass loss during the AGB phase introduces a second 
major source of uncertainty, and at the present moment only the comparison with observa- 
tions is expected to provide constraints on these points. 

Despite first successful models for the formation of the ^^C pocket have already been 
presented (Hollowell & Iben 1988; Herwig et al. 1997; Langer et al. 1999), the mass involved 
and the profile of ^^C formed at TDU in the He intershell are still to be assumed as free 
parameters, being related to hydrodynamical mixing mechanisms that cannot be properly 
accounted for in stellar models (Straniero et al. 2000). However a series of constraints can be 
obtained by comparing spectroscopical abundances in Pop. I AGB stars and in Pop. II Ba 
and CH stars (e.g.. Smith & Lambert 1990; Luck & Bond 1991; Vanture 1992; Plez, Smith, 
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& Lambert 1993; van Winckel & Reyniers 2000) with model predictions. In fact, Galactic 
disk AGB stars show a significant spread in s-process efficiencies that can be attributed 
to variations in the ^^C amount in the pocket and to the large sensitivity of the resulting 
s-process distribution on the initial metallicity. As discussed by Busso et al. (1995, 2000), 
the observed data are distributed around an average s-process efficiency rather similar to 
that of the main component in the solar system. However, the considerable intrinsic scatter 
revealed by spectroscopic studies can only fix broad limits to the amount of ^^C burnt. Notice 
that the same spread in ^^C concentrations was found to be appropriate to account for the 
s-process isotopic signatures of heavy elements in presolar grains recovered in meteoritic 
material (Zinner 1997), most likely condensed in the circumstellar envelopes of AGB stars of 
close to solar metallicity (Gallino et al. 1993, 1997; Lugaro et al. 1999; Amari et al. 2000). In 
our calculations the intrinsic spread in the s-process yields at each metallicity was modelled 
paramctrically by varying the ^^C concentration in the pocket by factors from to 2 times 
the "standard" (ST) value discussed in Gallino ct al. (1998) and Arlandini ct al. (1999) (~ 4 
10~^ Mq of ^^G). The ST choice was shown to reproduce the main component through AGB 
models in the mass range 1.5 — 3 Mq and metaUicity [Fe/H] = —0.3. Identical best-fits to 
the solar component are obtained for other choices of the ^^G pocket, provided the metaUicity 
varies according to the rule N(}^C) / Ni^^Fe) — const. Implications for different choices of 
the ST profile and/or metaUicity have been discussed elsewhere (Busso et al. 1995, 2000). 

Following a similar approach, Travagho et al. (1999) analyzed the GGE of Ba, La, Ge, 
Nd, Pr, Sm, and Eu. Both here and in Travagho et al. (1999) unweighted yields averaged 
over the whole set of ^^G abundances have been adopted. In Fig. 1 the abundances by 
mass of ^''^Pb in the He intershell material cumulatively mixed with the surface of a 1.5 
Mq star by TDU episodes are shown as a function of metallicity The main difference in 
the present computations with respect to those presented in Travagho ct al. (1999) lies in a 
finer interpolation grid at the lowest concentrations of ^^G. This larger grid of models can be 
useful for a better account of recent observational constraints from evolved stars at low and 
intermediate metallicity (Ryan et al. 2000a; Travagho, Gallino, & Busso 2000). Actually, the 
extension of the ^^G parameterization towards the lowest concentration of ^^G has negligible 
infiuence on the average unweighted yields (see § 3). 

The strong dependence of the s-process yields on stellar metallicity is evident in Fig. 1. 
The trend of s-process elements, and of ^°^Pb in particular, the most abundant among Pb 
isotopes, can be understood as follows. The build up of heavy s-process nuclei is due to 
neutron captures starting on pre-existing seeds; it is therefore formally of secondary origin 
and one would expect them to decline with declining metallicity. However, the abundance 
distribution is not only dependent on the initial Fe concentration, but also on the neutron 
addition efficiency, i.e. on the neutron exposure r (where r = / nnVrdt, with n„ neutron 



-6- 



density and vt thermal velocity) . The gradual increase of the neutron exposure toward low 
metallicities entirely masks this expected secondary behaviour (as first noticed by Clayton 
1988), resulting in a rather complex dependence of s-process yields on metallicity. Start- 
ing from AGB stars of nearly solar metallicity and going toward more metal-poor stars, 
the neutron-flux first builds up the s-elements belonging to the Zr-peak (at neutron magic 
number = 50). Then the Zr-peak is partly bypassed and the production of the elements 
at the second s-peak increases {N — 82, Ba-peak), reaching a maximum at [Fe/H] ~ —0.6 
(Travagho et al. 1999; Busso et al. 1999). For even lower metal contents, the n-flux feeds Pb 
(in particular ^'^^Pb), with a maximum production yield at [Fe/H] ~ —1 (see Fig. 1). The 
value of [Fe/H] at which the maximum is reached depends on the choice of the ^^C-pocket. 
For lower ■'■^C concentrations the production peak is shifted towards lower metallicities. Fi- 
nally, when very low metallicities are reached, also Pb starts to decrease, but it never really 
follows the expected secondary behavior. In fact, at very low metallicities, after the shortage 
of Fe during ^^C consumption becomes a conspicuous effect, the s-elements can continue 
to be fed thanks to the neutron captures starting on the abundant lighter primary neutron 
absorbers (e.g. ^^O, ^^Ne and its progeny). 

We also considered in the GCE model the role of intermediate mass AGB stars. In Fig. 2 
we show a plot similar to Fig. 1, but for the 5 Mq {upper panel) and the 7 Mq {lower panel) 
models. In these stars, the ^^Ne(Q!,n)^^Mg reaction is efficiently activated (Iben 1975; Truran 
& Iben 1977), since the temperature at the base of the convective pulse reaches values of 
T — 3.5 10^ K. As a consequence, the neutron exposure by the ^^Ne neutron source becomes 
more significant. Also the peak neutron density during the TP phase is consistently higher 
than in LMS {p ~ 10^^ n/cm^, see Vaglio et al. 1999, Straniero et al. 2000), overfeeding a 
few neutron-rich isotopes involved in important branchings along the s process, such as ^^Kr, 
®^Rb, and '^^Zr, with respect to the production of the s-only nuclei. As for the choice of the 
^^C neutron source, due to the much shorter interpulse phases in these stars (~6500 yr for a 
5 Mq and ~1500 yr for a 7 Mq) with respect to LMS (~ 3 — 6 10*^ yr), the He intershell mass 
involved is smaller by one order of magnitude. Consequently, also the TDU of s-process-rich 
material from the He intershell into the surface is reduced, again by roughly one order of 
magnitude. Due to the above reasons, for the 5 Mq and 7 Mq cases we considered as a 
new "standard" choice (STIMS) a ^^C mass scaled accordingly (M(^3C)stims = 10"'^ Mq). 
Then, in the calculations we considered a spread of ^^C abundances by a factor ±2 (from 
STIMS X 2 down to STIMS/2), as shown in Fig. 2. On the whole, IMS play a minor role 
in the production of Galactic Pb. A more extended set of evolutionary calculations are 
needed for IMS, however the model of 5 Mq and [Fe/H] = —1.3 does not show any relevant 
differences with respect to the one computed at solar metallicity (Straniero et al. 2000). 
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TABLE 2 

PB abundances in 10 METAL-POOR STARS 



Star 


[Fe/H] 


[Ba/Fe] 


[Pb/Fe] 


HD 


2665 


-1.95 


-0.32 


<1.0 


HD 


6229 


-1.01 


-0.04 


0.3: 


HD 


19445 


-1.97 


-0.08 




HD 


21581 


-1.57 


-0.06 


<0.5 


HD 


23439A 


-0.99 


+0.21 


0.6 


HD 


175305 


-1.31 


+0.01 


0.3: 


HD 


194598 


-1.05 


-0.02 




HD 


201891 


-1.04 


-0.04 


<0.5 


BD 


+23 3912 


-1.40 


+0.13 




BD 


+29 0366 


-0.92 


-0.08 


0.3: 



3. Lead abundances in stars 

Lead has very few lines suitable for an abundance analysis in the optical range. The 
only relatively unblended line in the solar spectrum is the Pb I line at 3683.48 A (Grcvcsse 
1969; Hauge & S0rli 1973; Youssef & Khalil 1989). A resolution R > 50,000 is required to 
adequately separate it from the nearby Fe I blend at 3683.623 A. Few spectra of high S/N 
for metal-poor stars have been taken in this difficult spectral region at this resolution. A 
line of similar strength is located at 4057.81 A; however this line is more heavily blended 
(mainly with CH lines), leading to more uncertain abundances. 

Hence, very few abundance determinations exist for Pb. Sncdcn et al. (1998) used UV 
spectra obtained with HST to derive Pb abundances in a few halo stars with large excesses of 
n-capture elements. The Pb I hue at 2833.05 A was clearly seen in HD 126238; a tentative 
detection was obtained in HD 115444, while the hne was not detected in the "normal" star 
HD 122563. Sneden et al. (2000) measured a Pb abundance from Keck spectra of CS 22892- 
52. Finally, Aoki et al. (2000) and Ryan et al. (2000a) recently gave a Pb abundance for the 
C-rich and s-rich star LP 625-44 from the line at 4057.81 A. 

To improve this situation, we used a set of spectra obtained in the last years at McDonald 
Observatory for another program (Gratton et al. 2000). These spectra have a resolution of 
R ~ 60,000. For ten of the stars considered in that program, the observed spectral range 
extended down to ~ 3678 A, allowing observation of the best Pb I line at 3683.48 A. These 
spectra are shown in Fig. 3, where they are compared with synthetic spectra obtained with 
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the atmospheric parameters listed by Gratton et al. (2000), and three different abundances 
of Pb: [Pb/Fe] = 0.0, 0.5, and 1.0. A mark signs the location of the Pb 1 line. In our spectra 
this is clearly out of the wings of the Fe I blend at 3683.623 A. Table 2 lists the abundances 
we inferred from these spectra; additionally, we also give abundances for Ba, derived from 
the same spectra. In most cases Pb abundances are uncertain, since the noise level is quite 
high at these wavelengths that are at the UV extreme of the spectra. Pb is clearly detected 
in the spectrum of HD 23439A, a star possibly classified as a CH star (Tomkin & Lambert 
1999), with small excesses of C and Ba ([C/Fe] = +0.09, [Ba/Fe] = +0.21). For this star 
we derive a quite large excess of Pb: [Pb/Fe] = +0.6. We think the Pb fine is also present 
in the spectra of HD 6229, HD 175305, and BD +29 366, although in all cases detection is 
uncertain. We did not detect the Pb line in the spectra of the remaining stars; upper limits 
were obtained for HD 2665, HD 21581, and HD 201891. 

Notice that according to Tomkin & Lambert (1999) HD 23439A belongs to a binary 
system of which HD 23439B (a single- lined spectroscopic binary) is also a dwarf CH star with 
a quite similar s-element composition. According to the authors the s-element enhancements 
must be primordial, reflecting a heterogeneous local interstellar medium. Unfortunately this 
companion does not belong to the sample of stars by Gratton et al. (2000). Analysis of the 
Pb I line in a spectrum in this star of comparable quality is highly desirable. 

4. Galactic chemical evolution of Pb 

The GCE model adopted here is described in detail in Ferrini et al. (1992) and Travaglio 
et al. (1999). As mentioned, it is based on the interconnected evolution of three zones: 
halo, thick disk and thin disk, whose relative composition in stars, gas phases, and stellar 
remnants is followed during the Galactic age. Here we consider the evolution of the solar 
annulus, located at 8.5 kpc from the Galactic center. The Star Formation Rate (SFR) is the 
outcome of self-regulating processes occurring in the molecular gas phase, either spontaneous 
or stimulated by the presence of massive stars. In Fig. 4 (small box) the SFR is plotted as 
a function of [Fe/H] during the evolution of the Galaxy (see also Travaglio et al. 1999). 

When introducing the s-process stellar yields of AGB stars only, we get the results listed 
in Table 3 (sixth column) for the Galactic s-process contribution to the various Pb isotopes 
at the epoch of solar system formation. In particular, att — tQ we obtain an s-fraction to Pb 
of 91%. The uncertainty in this evaluation may depend on the set of prescriptions adopted 
to estimate the s yields from AGB stars with varying the metallicity and on the general 
prescriptions adopted in the GCE model. However, as recalled before, a high degree of 
consistency is inferred from the detailed comparison of the predicted yields for the elements 
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TABLE 3 

S-PROCESS FRACTIONAL CONTRIBUTIONS AT t=t0 
WITH RESPECT TO SOLAR SYSTEM ABUNDANCES 





Solar(^) 






Strong 


(d) GCE(^) 




(%) 


(%) 


(%) 


(%) 




204pb 


1.94 


94 


79 


2 


90 


206 pjg 


19.20 


58 


31 


2 


60 


207pb 


20.62 


64 


30 


3 


77 


208pb 


58.30 


34 


10 


56 


89 


Pb 




46 


19 




91 



- Anders & Grevesse (1989) 

(b) _ Arlandini et al. (1999), "stellar component" 

(c) _ Arlandini et al. (1999), "classical component" 

- Beer, Corvi, & Mutti (1997) 
(e) _ xhis work 



from Ba to Eu with spectroscopic determinations in unevolved stars in the Galaxy (Travaglio 
et al. 1999). Concerning Pb itself, its s-fraction to solar depends on the present estimate 
of solar Pb based on meteoritic studies, which, according to Anders & Grevesse, is given at 
7.8% (one s.d.). In turn, the effect of the uncertainty of neutron capture cross sections on 
Pb isotopes is very small, because of the high precision of those determinations (see Bao et 
al. 2000). 

In the second column we also report, as reference, the elemental solar abundance of lead 
from Anders & Grevesse (1989). In the third column we show the prediction by Arlandini 
et al. (1999), as derived from the best-fit to the main component of the s-proccss in the 
solar system with the stellar model. The best-fit reported in this column was obtained with 
AGB models of metallicity [Fe/H] = —0.3 and the ST choice, by averaging the production 
in the mass range between 1.5 and 3 Mq. In the fourth column the updated predictions 
by the classical analysis are given (from Arlandini et al. 1999). Finally, in the fifth column 
we report the prediction from Beer, Corvi, & Mutti (1997) for the strong component of the 
phenomenological approach. 

From a comparison of the best fit by the classical analysis to the main component by 
Arlandini et al. (1999) with the outcome of our GCE, one can notice that there is a good 
agreement for all Pb isotopes, with the exception of ^'^^Pb. In fact Arlandini et al. (1999) 
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obtained an s-process contribution of 10% to ^°^Pb (19% to elemental Pb). They instead 

estimated a higher 34% of the s-process contribTition to ^'^^Pb (and 46% to elemental Pb) from 
their single AGB model ("stellar model") with the ST choice and [Fe/H]= —0.3 (actually 
a mean of 1.5 and 3 Mq AGB models), which also reproduces the main component. The 
discrepancy of the two above predictions on ^'^^Pb by the classical analysis and by the "stellar 
model" was already a clear indication that not a unique AGB stellar model, nor the classical 
analysis, were able to explain the main component in the solar system, which must be 
considered as the outcome of different generations of AGB stars prior of the solar system 
formation. 

As a matter of fact, our GGE calculations provide 89% of ^osp]^ abundance, thanks to 
the contribution of different generations of AGB stars and in particular to those at low metal- 
licities, which are the main contributors to ^°^Pb. The proper account of GGE confirms in a 
quantitative way what was anticipated by Gallino et al. (1998), i.e. that the role previously 
attributed to the strong component is actually played by AGB stars of low metallicity. 

A remarkable result of introducing the AGB yields at different metaUicities (as shown in 
Fig. 1 for ^°^Pb) as input data of GGE computations is that the predicted chemical compo- 
sitions in the interstellar medium at any epoch are grossly independent of the particular ^^G 
profile adopted inside the ^^G pocket. For example, the results presented in the last column 
of Table 3 remain essentially the same even assuming a constant ^^G mass fraction in the 
pocket. A more detailed analysis will be presented elsewhere, including the GGE trends of 
all Pb isotopes. 

The s-process yields discussed in § 2 and shown in Fig. 1 allow us to estimate the 
chemical enrichment of Pb in the Galaxy at different epochs. We show in Fig. 4 the resulting 
s-fraction of lead in the thick and in the thin disk. This has been obtained by considering 
AGB stars in the mass range 2 — 8 Mr^, but the dominant production of Pb comes in the 
mass range 2 — 4 Mq (see discussion in § 2). It is evident that the s-process contribution 
dominates the GGE of Pb, starting from [Fe/H] ~ —1.5, both in the thick and thin disk. 
This is mostly due to the coincidence between the peak of Pb production in LMS (in the 
range —1.5 < [Fe/H] < —1) with the peak of the star formation rate in the thick and thin 
disk. At lower values of [Fe/H], independently of the characteristics of the GGE model, the 
contribution of s-process nucleosynthesis rapidly decreases as a result of the long lifetimes 
of AGB stars. 

This confirms previous results of the GGE evolution of elements from Ba to Eu (Travaglio 

et al. 1999). At very low mctallicitics, as suggested by spectroscopic observations (in partic- 
ular of Ba and Eu) and pointed out by Travaglio et al. (1999), the s-contribution by AGB 
stars is by far too small. Instead, as anticipated by Spite & Spite (1978) on observational 
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grounds and by Truran (1981) on theoretical grounds, the heavy element abundance pattern 
in very metal-poor stars is compatible with an r-process origin. This point was recently 
sustained by observations of low metallicity halo stars (Sneden et al. 1998; McWilliam 1998; 
Sneden et al. 2000; Westin et al. 2000; Ryan et al. 2000b, and references therein). 

Prom the theoretical point of view, despite a large number of recent works, the astrophys- 
ical site of the r-process is still a debated problem (e.g.. Wheeler, Cowan, & Hillebrandt 1998; 
Frciburghaus et al. 1999). In order to quantify the r-contribution, as described in more detail 
by Travaglio et al. (1999), we treat the r-process as a typical primary mechanism occurring 
in a subset of Type 11 SNe (those in the mass range 8 — 10 Mq). Recent theoretical simu- 
lations by Wheeler et al. (1998) support this choice. Our estimate of r-process abundances 
at t ^ ^0 was simply derived by subtracting from the solar abundances the s-fractions. 

In Fig. 4 the predicted Galactic evolutionary trend for the s + r contribution to Pb is 
shown. In the figure we also compare our model to recent spectroscopic data for halo and 
disk field stars at different metallicities, discussed in the previous Section. It is clear that 
a larger set of observations is desirable, and we underline that these are the first available 
observational data of Pb. The star with [Fe/H] = —3.1 and [Pb/Fe] = -1-1.3 is the very 
pecuhar CS 22893-052 (recent observations of Pb in this star has been presented by Sneden 
et al. 2000). This star shows exceptionally high abundances of all heavy elements, with 
r-process/Fe enrichment of ~40 times the solar ratios. Therefore it was probably born in an 
environment strongly polluted by SN debris with respect to the average interstellar gas in the 
halo. This supports the current idea of an incomplete mixing of the gas in the Galactic halo at 
early epochs, which would allow the formation of stars both very deprived, or very enriched, 
in r-process elements (for recent studies of this problem see e.g. Tsujimoto, Shigeyama, & 
Yoshn 1999; Raiteri et al. 1999; Argast et al. 2000; Travagho, GaUi, & Burkert 2001). 



5. Conclusions 

In this paper we have studied the evolution of Pb in the interstellar gas of the Galaxy, 
adopting a new set of models for s-processing in AGB stars of different metallicities, and 
new observational constraints of unevolved field stars for Pb abundances. 

We have calculated the s-process yields with post-process calculations based on AGB 
models and applied these yields in the framework of a Galactic model that follows the 
GCE of the halo, thick and thin disk. Taking into account the important role played by 
different generations of low-metaUicity AGB stars for the production of Pb in the Galaxy, 
we obtain at i = a Pb s-fraction of 91%, together with a sharp increase at [Fe/H] ~ —1.5. 
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Confirming the prediction by GaUino et al. (1998), we also found that the production of the 

most abundant Pb isotope, ^°®Pb, is matched when the s-process occurring in low-metaUicity 
AGB stars is properly considered in the context of the chemical evolution of the Galaxy. Only 
a small portion of ^^^Ph can be attributed to the r-process, even accounting for decays of 
transuranic isotopes. 

We are indebted to F. Ferrini for allowing us free access to his GCE model, to O. 
Straniero, A. Chieffi and M. Limongi for sharing the results of FRANEC computations, and 
for a longstanding collaboration on scrutinizing deeper and deeper in AGB stellar evolution- 
ary models. We also thanks S.G. Ryan for very useful discussions. 
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Fig. 1. — Abundances by mass (relative to solar) of ^°^Pb in the He intershell material 
cumulatively mixed with the surface of a 1.5 Mq star by third dredge up episodes as function 
of metallicity, for different assumptions on the mass of the ^^C pocket. The thick continuous 
line represents the unweighted average of all cases shown. 
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Fig. 2. — The same of Fig. 1, but for a 5 Mq star {upper panel), and for a 7 Mq star {lower 
panel). In both panels the thick continuous lines represent the unweighted average of the 
cases shown. 
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3683 3683.5 3684 3683 3683.5 3684 3684,5 

Wavelength (A) Wavelength (A) 

Fig. 3. — Comparisons between observed spectra for ten metal-poor stars (open squares) and 
synthesized spectra for the same stars {solid lines) for a small spectral region including the 
Pb I line at 3683.48 A. Synthetic spectra were computed with the atmospheric parameters 
listed by Gratton et al. (2000), and three different abundances of Pb: [Pb/Fe]=0.0, 0.5, and 
1.0. A mark signs the location of the Pb I line. 
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Fig. 4. — Galactic evolution of [Pb/Fe] according to our Galactic model predictions, including 
both s— and r— process contribution for the halo {short- dashed line), thick disk {long-dashed 
line), and thin disk {solid line). The Pb s-fraction from AGB stars in the thin disk is also 
shown with thick-dotted line. Observational data are from Sneden et al. (1998) {filled circle), 
Sneden et al. (2000) {filled triangle), and from the present work {filled squares). Error bars 
are shown only when reported for single objects. In the lower-right small box is also shown 
Star Formation Rate (in unit of Mq pc~^ Gyr~^ and logarithmic scale) for halo, thick disk, 
and thin disk, according to our model predictions, as a function of [Fe/H]. 



